Is Eddington-Born-Infeld theory really free of cosmological singularities? 
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The Eddington-inspired-Born-Infeld (EiBI) theory has been recently resurrected. Such a theory is 
characterized by being equivalent to Einstein theory in vacuum but differing from it in the presence 
of matter. One of the virtues of the theory is to avoid the Big Bang singularity for a radiation filled 
universe. In this paper, we analyze singularity avoidance in this kind of model. More precisely, we 
analyze the behavior of a homogeneous and isotropic universe filled with phantom energy in addition 
to the dark and baryonic matter. Unlike the Big Bang singularity that can be avoided in this kind 
of model through a bounce or a loitering effect on the physical metric, we find that the Big Rip 
singularity is unavoidable in the EiBI phantom model even though it can be postponed towards a 
slightly further future cosmic time as compared with the same singularity in other models based on 
the standard general relativity and with the same matter content described above. 



Einstein theory of general relativity (GR) is an ex- 
tremely successful theory for nearly a century [l[. De- 
spite of all its advantages, it is expected to break down 
at some point at very high energies, for example in the 
past evolution of the Universe where the theory predicts 
a Big Bang singularity and the laws of physics cease 
to be valid. This is one of the motivations for looking for 
possible extension of GR. In addition, it is hoped that 
modified theories of general relativity, while preserving 
the great achievements of GR, would shed some light 
over the unknown fundamental nature of dark energy or 
whatsoever stuff that drives the present accelerating ex- 
pansion of the Universe (see for example Ref. 0] ) . 

There have been many proposals for alternative theo- 
ries of GR almost as old as the theory itself. One of the 
oldest was proposed by Eddington where the con- 
nection rather than the metric plays the role of the fun- 
damental field of the theory. The gravitational action 
proposed by Eddington back in 1924 [3| is equivalent to 
Einstein theory of GR in vacuum. One of the weak points 
of the theory is that it does not incorporate matter. Re- 
cently, an Eddington-inspired-Born-Infeld theory (EiBI) 
has been proposed in Ref. [H where matter fields are 
incorporated into the Lagrangian formulation. More im- 
portantly, it turns out that this theory avoids the Big 
Bang singularity that would face a radiation-dominated 
universe in standard GR jsj. The apparent fulfillment 
of the energy conditions in EiBI theory was considered 
in Ref. @, where the adjective apparent refers to quan- 
tities defined with respect to a metric compatible with 
the connection that defines the theory. Their analysis 
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leads to a sufficient condition for singularity avoidance. 
Besides, the gravitational collapse of nonintcracting par- 
ticles, i.e., dust or equivalently pressureless matter, does 
not lead to singular states in the nonrclativistic limit 
(Newtonian regime) Q (see also Q). This theory has also 
been studied as an alternative scenario to the infiationary 
paradigm Q. Furthermore, possible constraints on the 
parameter characterizing the theory have been obtained 
using solar models [loj . neutron stars [ll| . and nuclear 
physics [13]. It has also been shown that such avoidance 
of Big Bang singularity is more general and not limited 
to the radiation-dominated universe [l3|. Despite of all 
the virtues of the EiBI theory, a cosmological tensor in- 
stability in this model was found in Ref. [ij] . In addition, 
this theory behaves similarly to the Palatini f(R) grav- 
ity and shares the same pathologies, such as curvature 
singularities at the surface of polytropic stars and some 
unacceptable phenomenology |la |. 

In this letter, we ask the simple questions: Is EiBI 
theory really free of cosmological singularities? In par- 
ticular, is the theory free of dark energy related singu- 
larities? In Ref. @, it was shown that if the null en- 
ergy condition is fulfilled, then the apparent null energy 
condition is satisfied. It turns out that the null energy 
conditions are not always fulfilled; a clear example of 
it is a super-inflationary phase within GR. Moreover, in 
recent years a new singularity named the Big Rip has 
been identified where the null energy condition is in fact 
not fulfilled and the Universe is ripped apart: the Hubble 
rate and its cosmic derivative approach infinity in a finite 
cosmic time [l^ [l3l . Can such singularity be avoidable 
in the theory proposed in Ref. @ ? This question is even 
more pertinent in the aftermath of the release of WMAP9 
data, which hints on the possibility of a phantom energy 
component in the Universe more pronouncedly than that 
deduced from the WMAP7 data [l^. The analysis of 
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the possible occurrence of a Big Rip in the future of the 
Universe is therefore timely. 

Our starting point is the gravitational action with the 
metric g^j/ and connection recently proposed in Q : 



5n.(.g,r,M') 



(1) 



where R^,^ (F) stands for the symmetric part of the Ricci 
tensor and, as indicated in Eq. ([1]), is constructed from 
the connection F. We consider the action under the Pala- 
tini formalism, i.e., the connection F^^^ is not the Levi- 
Civita connection of the metric 5^^. The parameter n is 
a constant with inverse dimensions to that of the cosmo- 
logical constant (in this letter, we will work with Planck 
units SttG = 1 and set the speed of light to c = 1), A is 
a dimensionless constant and 5m(5,F, stands for the 
matter Lagrangian. This Lagrangian has two well defined 
limits: (i) when is very large, we recover Eddington's 
theory and (ii) when \kR\ is small, we obtain the Hilbcrt- 
Einstein action with an effective cosmological constant 
A = (A — 1)/k d]. A solution of the action in Eq.(IT|) can 
be characterized by two different Ricci tensors: Rf^^{T) 
as presented in Eq.ljT]) and R^j,u{g) constructed from the 
metric g. There are in addition three ways of defining 
the scalar curvature. These are: g^^''R^y{g), g^^'^ R^yiT) 
and i?(r). The third one is derived from the contraction 
between i?^i/(r) and the metric compatible with the con- 
nection F. Therefore whenever one refers to singularity 
avoidance, one must specify the scalar curvature(s). 

The equations of motion are obtained by varying the 
action, ^EiBij with respect to the metric and the con- 
nection. The energy-momentum tensor is conserved in 
this theory. Thus, for a Friedmann-Lemaitre-Robertson- 
Walker (FLRW) universe filled with a perfect fluid with 
energy density p and pressure p, we obtain the familiar 
relation 



p + m{p + p) = 0. 



(2) 



After the variation of the action ([T]) and combining it 
with the conservation equation Eq.Q, we arrive at a 
modified Friedmann equation for a universe with scale 
factor a, which is filled with a perfect fluid with energy 
density p and pressure p Q : 



where pt 



A, Pt =p- A, A = (A- 



F = 2- 
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U ^{1- Kpt)Hl + Kpty^^V^ + Kpt)^K 
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(4) 
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Consequently, for a universe whose matter content is 
dominated by a single component with its equation of 



state (EOS) p = wp, and in the absence of a cosmologi- 
cal constant, i.e., A = 0, its evolution is governed by 



3k 



(1 + 3w)p -2 + 2^/(1 + P)(l 
(4 + (1 - w){l - 3w)p + 2w{l + 3w)p2)2 



(6) 



where p = up [l3|. It can be easily verified that the 
Big Bang singularity can be avoided in this theory for a 
radiation dominant universe; i.e., w = 1/3. More specifi- 
cally, the Universe either bounces in the past for the case 
of K < 0, or has a loitering behavior in the infinite past 
for the case of k > Q. Despite that the Big Bang 
singularity is avoided with respect to the metric g; i.e., 
the Hubble rate, its cosmic time derivative, scalar curva- 
ture g^'^ R^y{g), and Ricci curvature R^y{g) are all finite, 
nevertheless the scalar curvature of the metric compati- 
ble with the connection F, i.e., i?(F), still diverges when 
the scale factor approaches the minimum (see TABLE 

A natural question, inspired by this finding, is how 
general is the singularity avoidance in EiBI as compared 
with GR. In particular, can EiBI cure or smoothen the 
Big Rip singularity? Such a singularity is expected in 
GR for a phantom energy dominated universe with a con- 
stant equation of state. In order to address this question, 
we focus on the late-time evolution of a FLRW universe 
filled with phantom energy (w < — 1 and is constant) in 
addition to the dark and baryonic matter. The matter 
content reads 



Pt 
Pt 



Pm + Pi 

Pw = wpw^a 



Pn 



-3(l+w) 



-3(1+M>) 



(7) 



where p^ and pw are the density of matter and dark 
energy, respectively. 

This model contains four parameters, Pmoi Piooj ^- ^^'^ 
K, but only three arc independent because of the cosmo- 
logical constraint obtained by evaluating the Friedmann 
equation at the present time, which read^: 



where 



2W 



(8) 



X = 2- 
Y 



2AB ' 

n^[i + w{i + w)fi^{fi^ + fi^)] 
c 



(9) 



and A = 1 - wHk^ui, B = I + il^{ilm + ^w), C 
^rn + (1 + w)r2tu, respectively. 



^ The parameters Q,rn and are defined in the standard way, i.e., 
f2m = pmal Pc and Cm = /Pc, where pc is the critical density. 
In addition, VLk, = ftpc = ZHq^k, whore Hq is the current Hubble 
parameter. 
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At small r^K, f{^K,^m,^w,'U}) can be written as 
O^, n^,w) K {iljn + ^w)^H, + 0(r2«^), which con- 
firms again that EiBI reduces to GR for vanishing £7^. 

Since EiBI theory contains a new parameter, k, as com- 
pared with GR, we shall first put some constraints in it 
before proceeding further. The estimation will be based 
on three points as follows: (i) We expand the right-hand- 
side of the cosmological constraint, Eq.®, up to the sec- 
ond order in and obtain the solutions: 



simply combining Eq. ([2]) and Eq. (|3]): 



f^K =0, 

oiv 



(10) 



where 



(11) 



This approximation is used here as a way to simplify 
the presentation of our results and does not affect the 
conclusions of our paper, (ii) We assume that the model 
conforms with the u^CDM scenario at present so as to be 
able to explain the current acceleration of the Universe. 
We may therefore assume that 0.267 < flm < 0.287, 

0. 713 < flyj < 0.733, and -1.147 < «; < -1.021 [H. 
These cosmological constraints imply that K < (cf. 
Eq. ([TT]) ). (iii) We restrict our model to a positive r2„, 

1. e., K > 0, in order to avoid the imaginary effective sound 
speed instabilities usually present in EiBI theory [Il| . As 
a consequence, fl^ vanishes whenever rim + ^w < 1- (See 
Eq. ^ and (ii).) 

In summary, under the above conditions we can con- 
clude the following: (i) 17^ = whenever -I- fiu, < 1, 
where we recover GR. (ii) For ftm + > 1, we may 
consider the second solution in Eq. (fTO)) . which is posi- 
tive in this case. In the latter case, we will assume that 
is small, i.e., the deviation of EiBI theory from GR 
is small, so that Vim + ^ 1 is iu agreement with the 
observational data [l^ . (iii) One can always find a suit- 
able value for or k, to fit a specific set of parameters 
firm and w. 

We now investigate the asymptotic behavior of the 
Universe within this framework. This amounts to de- 
termining the Hubble parameter, H , and its cosmic time 
derivative, H, at large scale factors. From here on we 
set w = —1 — e, where e is positive. As the dark energy 
corresponds to the phantom matter in our setup and the 
Universe is expanding, the conservation of such an energy 
density implies a growth of p^, (see Eq. ([2])), unlike the 
baryonic and dark matter, which would quickly become 
negligible as compared with p,„ . We therefore neglect pm 
in our estimation of H and H . 

Under the above assumptions, (1 <C Kpt ~ np^), we 
obtain the asymptotic behavior of the Hubble parameter 
H given in Eq. ([3]) and the cosmic time derivative, H, by 



4v^ 



3(2 -h3e)- 



-Pt, 



(2 + 3.)2 'P'- 



(12) 



The above results correspond to the dominant terms in 
the expansion of kH^ and kH as functions of npt ~ upw 
As can be seen, H and H will blow up when np^ diverges 
at an infinite radius of the Universe. 

We can also prove that a phantom energy dominated 
EiBI universe has a well defined for any value of p^ . 
In fact, the square of the Hubble parameter in Eq. (|6]), for 
p — Kpw, is positive-definite and it vanishes only when 
Pw = 0. 

So far we have shown that the total density pt, the 
total pressure pt, H^, and H will all diverge when the 
scale factor goes to infinity. Our next step is to confirm 
the existence of the Big Rip singularity at some finite 
cosmic time. The cosmic time can be evaluated directly 
from the integral 



i^o(ising — to) 



dz 



_i {l + z)Eiz)' 



(13) 



where z is the redshift parameter, E{z) — H/Hq, iging 
and tg are the cosmic time the singularity takes place 
and the present time, respectively. For the relativistic 
Friedmann equation, we have Ec\{z) = ['[lm{z + 1)^ + 
ilw{z + and for the EiBI modified theory, we 

can use Eq. (|3]) with the matter content given in Eq. ([7]) 
to derive EEmi{z), where E^iBiiz) = i?EiBi(z)/^^o- We 
show in TABLE [III the results of our numerical integra- 
tion based on Mathematica 7, where we have assumed 
flm = 0.287 and fl^ = 0.733 [l|], and used the con- 
straint Eq.® to find the corresponding $7^- We choose 
those limiting values of flm, to enhance the possible 
effects of the EiBI model; i.e., we choose observational 
values that maximize the inequality fim + ilw ^ 1. We 
see that the cosmic time derived from the EiBI theory 
is finite and of the value of ten times the current age of 
the Universe, which implies that this theory is not able 
to remove the Big Rip singularity occurring in GR, even 
though this singularity can be slightly pushed towards a 
future time as compared with GR. 

Our results indicate that the scalar curvature con- 
structed from the physical metric (/^^ will blow up at the 
Big Rip. It can be shown that i?^i/(r) and g'^''R^i,{T) 
also blow up at tging where a diverges, whereas R{T) re- 
mains finite. Specifically, i?oo(r) = (1 ^ U)/k — > — oo 
and R^J{^) = [a'^{V - l)Sij]/K oo, while i?(r) = 
(U - 1)/Uk + 3iV - 1)/Vk 4/n, as a oo. 

An interesting model for a modified theory of gravity 
was suggested in Ref. Q . It was shown that in this model 
the Big Bang singularity for a radiation-filled universe 
can be removed [5| , but the scalar curvature constructed 
from the metric compatible with F still blows up as we 
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have shown. On the other hand, it is known that for 
the class of dark energy models with w < —1, i.e., the 
phantom models, the Big Rip singularity is inevitable for 
a constant w in the framework of GR. Our main objec- 
tive of this paper is to see if the Bafiados-Ferreira EiBI 
model can help also to remove the Big Rip singularity. 
We tackled this issue by investigating the possible oc- 
currence or avoidance of doomsdays in this model. We 
analyzed an EiBI FLRW universe filled with dark matter 
and phantom energy with a constant equation of state. It 
is well known that a universe with such a matter-energy 
content under GR would face a Big Rip. Our result indi- 
cates that the Big Rip singularity remains inevitable in 
the EiBI theory albeit providing a minor postponement, 
as shown in TABLE |lll The onset of the Big Rip is in- 
dependent of the amount of dark matter or dark energy; 
i.e. fim and n^^. In fact, the scale factor, the Hubble 
parameter and its cosmic time derivative all blow up in 
a finite cosmic time. Consequently, the scalar curvature 
constructed from the physical metric g^i^ will also blow 
up. We have shown as well that i?^,y(r) given in the 
action in Eq.([T]) and g^'^R^y{T) are infinite at the singu- 
larity, whereas R{T) remains finite. The key message to 
take home from this letter is that a Big Rip singularity 
cannot be avoided in the EiBI model but it is smoother 
than that in GR. This is unlike the Big Loitering in a 
radiation dominant EiBI universe, which is rougher than 
that in GR, as shown in TABLE H] 

We will present elsewhere the behavior of other dark 
energy related singularities/events [l9l - [2]| such as big 
freeze, sudden singularity, type-IV singularity, little rip, 
etc., in the EiBI framework 123. 
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TABLE I. The behavior of different possible Ricci tensors and 
scalar curvatures in the EiBI model at the Big Loitering in a 
radiation dominated universe and at the Big Rip singularity 
in a phantom dominated universe. 



0.021 
0.041 
0.061 
0.081 
0.101 
0.121 
0.147 



Jyo(tsing-to)(GR) 



37.0149 
18.9291 
12.7039 
9.55371 
7.65167 
6.37884 
5.24246 



Ho(tsing -to)(EiBI) 



37.1153 
19.0294 
12.8041 
9.65379 
7.75166 
6.47876 
5.34228 
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TABLE II. The cosmic time elapsed from the present time 
to the Big Rip singularity time, normalized to the current 
Hubble parameter; i.e., -ffo(tsing — to), for different values of 
e in GR and in the EiBI theory. We see that such cosmic 
time remains finite in the EiBI theory, meaning that the Big 
Rip singularity is inevitable. Here we assume Q,m = 0.287 
and Q.-U, — 0.733, then use the constraint Eq. to find the 
corresponding Q^.. 
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